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FOREWORD 


This  report  was  prepared  for  Applied  Mathematics  Research  Laboratory, 
Aerospace  Research  Laboratories,  by  A.K.  Chattopadhyay  under  Project  7071, 
Research  In  Applied  Mathematics.  This  work  was  performed  at  U.S.A.F. 
Aerospace  Research  Laboratories  by  the  author  while  In  the  capacity  of 
Technology  Incorporated  Visiting  Research  Associate  under  contract 
F33615-71-C-1463. 

In  this  report  the  author  studies  the  robustness  of  Hotelling's 
2 

Generalized  Tq  test  under  the  violation  of  general  linear  hypothesis 
both  In  respect  of  mean  and  variance  covariance  matrices. 

The  author  wishes  to  thank  Dr.  P.R.  Krlshnalah  for  some  useful  dis¬ 
cussions,  and  Mrs.  Georgene  Graves  for  typing  the  manuscript  carefully. 


ABSTRACT 


In  this  paper  the  author  derived  asymptotic  expressions  for  percentile 
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and  c.d.f.  of  i lotel ling's  Generalized  T0  statistic  under  the  nonnull 
assumption  of  mean  matrix  and  variance  covariance  matrix  satisfy  (3)  and 
(4)  given  in  the  text. 

These  expressions  can  be  used  to  study  the  robustness  of  the  cest  with 
respect  to  power  function  and  stabilization  of  critical  region. 
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1.  Introduction 

In  an  earlier  paper  [2]  the  asymptotic  formulae  for  the  distribution 
and  percentile  of  statistic  T  ■  mtrS^S2‘^  have  been  obtained  up  to  terms 
^  where  mS^  and  ..$2  are  independently  distributed  W(m,p,E-,)  and  W(n,p,E2) 
respectively.  Similar  expansions  for  the  ratio  of  two  independent  trace 
statistics  are  also  obtained.  This,  in  fact,  generalizes  the  work  of 
previous  authors  [3],  [5],  [6],  In  a  recent  paper  [6]  Slotanl  gave  an 
asymptotic  expansion  for  the  nonnull  distribution  of  Hotelling's  general- 

p  *1 

Ized  T  up  to  terms  -U-  by  using  James  [4]  and  Welch  [7]  Idea  by  expanding 
0  rr 

the  characteristic  function  by  the  perturbation  technique. 

In  this  article  we  generalize  the  earlier  results  and  find  asymptotic 

expansion  up  to  terms  of  order  ^  for  c.d.f.  and  percentile  of  the  trace 

statistic  when  mS1  has  W(m,p,E,n)  and  E-*  *  Eg  but  otherwise  satisfy  (3)  and 

2 

(4).  The  expression  (6)  can  be  used  to  compute  the  power  of  TQ  test  for 
the  generalized  linear  hypothesis  when  departure  from  hypothesis  M  3  0  and 
Ei  ■  e2  Is  present. 

It  can  also  be  used  In  cases  to  test  two  covariance  matrices  when  one 
has  a  noncentral  Wlshart  distribution. 


s 


2.  Formulation  of  the  Problem 

Let  Z  =  {2^  ....  z^}  be  a  p  x  I"  random  matrix  where  z^'s  are  Indepen¬ 
dently  distributed  according  to  p  variate  normal  distribution  with  mean 
vector  vi.  and  variance  covariance  matrix  £,  =  B'^.  Let  nS  =  n(s -  .)  be  a 

1  in  ij 

p  x  P  matrix  which  is  independent  of  Z  and  follows  a  central  Wishart 
distribution  W(n,p,A"^)  with  n  degrees  of  freedom  and  variance  covariance 
matrix  -  A"  .  Hotelling’s  generalized  is  given  as 

7  -i 

V  =  tr  Sn  ZZ 

o  n 


m 


-1. 


■^^l  Sn  Zi  •  when  S1  *  E2 


2 

Our  aim  is  to  find  asymptotic  expansion  for  percentile  and  c.d.f.  of  TQ 
when  f  z2  but  otherwise  satisfying  (3). 
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3.  Asymptotic  Expansion  for  Percentile 
Lc  c 

r«(e)  =  Pr  {tr  Sn_1 ZZ  '<  0} 

Now 

Pr  { tr3Z7.  '  <  b) 

2 
OU 

=  e"?  l  [7  I  Ix^J-'eMx 


2' 

j 

w 

r 

.u 

j!^  er+T  J. 

Where 


*  x2mp  {0,a,2) 


w2  ■  trBMM*  =  trn 


M  *  (u-j ,  . . um>  f1  0,  p  *  m£ 


2  2 

and  x  mD  (0,  u  }  Is  the  c.d.f.  of  noncentral  chi-square  variable  with 

2 

mp  d.f.  and  noncentrality  parameter  w  . 

Now  we  can  try  to  find  a  function  h  (S_)  of  the  elements  of  S  n 

"  n  * 

belnq  large  enough  such  that 

(0)  -  Pr  {tr  Sn'1ZZ  <  h  (Sn)J 
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NOW 


G(o)  *  E  {Pr  {exp{tr($  -  A'1)  3} 
n  n 

Pr  {tr  AZZ'  <  h(A"*‘)})} 

■  »  Pr  {tr  AZz'  <  h(A-1)} 

Where 

.  _n 

0*  exp  [-tr  A'1  3]  |  i  .  1  a'1  3  |  ^ 

'  '  +  "  rLVrS  0t“  *st  8“r  *  0(n'i’)  H) 

HP  «  P)  ■  <%  0  *  «(J> 

°n  °ip 

•  • 

A'1-  . 

>1  ’  ’  ’  °PP. 

Now  expanding  h(Sn)  around  e  we  get 

h(Sn)  ■  e  ♦  hj{Sn)  ♦  h2(Sn)  +  ... 

Where  hi(Sfl)  Is  Ofn'1) 

Thus  expanding  h(Sn)  around  h(A~^)  we  qet 

GM  •  C  !st  V*  0("'2>I 

4 
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[1  +  h^A-1)  D  +  0(n"2)]  Pr  {tr  AZZ*  <  e  } 

Where  D  8 

o  0 

This  being  true  for  all  large  n  we  get 

[h,  (A'1)  D  ♦  i  I  °tU  3st  >J  Pr  (tr  AZz' <  6)  -  0 
Again  let 

J  =  Pr  {tr  (A-1  +  e)'}Zl  '  <  6} 

Following  [2],  [3],  [6]  we  get 

in  2 

J  =  |I  -  xa|"*  exp(-|  ) 

exp  {%  E  tr(I  -  xa)"1  fl}  X2mp(e»  0)  (2) 

Where 

a  *  E  -  1,  £r  x2mp(  o »  <*>2) 

“  x2^  +  2r  (e,  u.2)  and 

x  -  +  e)'1  -  I 

-  (B^A  -  I  )-[EfS  (8_1A)  (A’^A) 

+  £  crs  etu  (A_1rsA^  ^tu  A>  '  - 

Where  A"^  is  the  p  x  P  matrix  with  (1,j)th  element 


*<6r1  6si  +  6rj  5s1> 
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wm 


flow  let 

| chi  (B~ 1 A  -  I)|  *  | chi  (F) j  <1,1=1,  ...  p  (3) 

Where 

B-1A  -  I  *  F 

Expanding  (2)  and  equating  coefficients  of  ers  Et(J  with  those  in 
Taylor's  expansion  of  0  around  e  ■  0  and  denotlnn 

tr  (  A ~ 1  rsA )  **  (rs) 

tr  (A_1rsA)  (A-1tuA)  -  (rs|tu) 

tr  FtA'^A)  (A~^trA)  -  (F|rs|  tu) 

tr  F2  -  (F|F) 

tr  F  «  (F)  etc. 

and  using  the  following 

l  o$t  our  (rs  ] tu)  ■  Jjp(p  +  1) 

[  (rs)  »  p 

l  ost  our  (rs)  (tu)  -  p 

I  °st  aur  (Flrs*  *tu)  *  *F) 
l  ost  our  ( F | rs | tu )  •  (F)  (p  +  1)  /2 


I  ost  our  ( F ) rs )  ( F 1 1 u)  =  (F|F#) 
l  ost  our  (n I  rs  I F I  tu)  *  ’;[(n)  (F)  +  ( nF ) ]  etc. 
get  after  neglecting  terms  involving  f  when  F  *  (f^) 

-  h1  (A"1)  D  Pr  {tr  Ml '  <  0} 

1  4  2 

■  in  .L  ,j(m,p)  V  ♦  "  ) 

j  * 

+  In  j l0  bj("’p)  *W  +  2j(0>  “2) 

+  0(n'2) 

Where 

a0(m,p)  =  mp(m-p-l) 

2 

a-|(m,p)  =  -2m{mp-w  ) 

a2(m,p)  58  mp(m+p+l)  -  2(2m+p+l)w2  +  trn2 

a.j(m,p)  =  2{(m+p+l)w2  -  tm*"} 

2 

a4(m,p)  s  trn 
_2 

bq(m,p)  =  -  -j—  (F)  n(m-p-l) 

b^m.p)  -  5  (F)  I4”1  -  _  mp2  ‘  nn  "  2(n+n+1) 

-  ^(nF)  (n  (m-p-1)  +  4} 


(4) 
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/c\  22  2  2  3 

b2(m,p)  *  *  P  +  8m  ♦  m  p  +  3m  p  +  4mp  +  4m 


-  (6nr  +  8mp  +  8m  +  4)  (n)} 


-if-  {mp  +  mp2  -  3m2p  -  16m  -  4p  -  8  -  (4m  +  2p  +  2)  (n)} 


♦  2m(oF  )  -  2(n  Fn) 


b.j(mfp)  *  {m2p2  +  m2p  ♦  4mp  +  4m  +  m3p  +  4m2 


-  (6m  +  lOmp  +  10m  +  8)  (n)} 


-  (3m2p  +  mp  +  mp2  +  2m  +  12p  +  20 


-  (12m  +  8p  +  8)  (fi)>  ♦  6(n  Fp)  +  2(nFn  ) 


b4(m,p)  ■  {m2p  +  mp2  +  mp  +  (2m2  +  4mp  ♦  4m  +  4)  (n)} 
♦  Jga.  (8m  +  8p  ♦  12  -  (lOp  +  22)  (n)} 


+  2m(nF  )  -  6(n  Fn) 


b5(m,p)  »  2(nF)  (n)  (m  +  P  ♦  1)  +  2(n  Fn) 


+  2(nFn  ) 
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Where  as  noted  earlier  we  dropped  terms  Involving  f  f.  ,  etc.,  and 

2  j  k  I 

9mp{0,w  )  means  the  c.d.f.  of  noncentral  chi-square  distribution  with  mp 

2 

d.f.  and  noncentral Ity  oarameter  «  . 

Thus 

h<sn>  ’  •-  V2j<5,“2) 

*Tn  j„  bj<m*P)  V2J(a‘“2)3  [G' 

*  0(n’2)  (5) 

Where  e  Is  the  corresponding  percentile  of  linear  function  of  non¬ 
central  chi-square  variable  of  the  form  [1] 

v  ■  j,  y/(m) 

Where  x j ' s  are  the  characteristic  roots  of  AB~^  and  G(e)  Is  the  c.d.f. 
of  Y. 

Here  our  form  (5)  differs  slightly  from  those  given  In  [2],  [3],  [5], 
[6]  but  this  form  gives  a  uniform  result  both  for  the  percentile  and  for 
the  c.d.f.  as  given  below. 
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4.  Approximation  for  c.d.f.  of  Hotelling's  Generalized  T2  Under  Nonnull 
Assumptions  on  Mean  Vector  and  Variance  Covariance  Matrix 

Here  we  proceed  as  in  [2],  [3],  [5]  and  using  our  earlier  calculation 
we  get 

Pr  {tr  S^1  11  <_  0}  *(H)Pr  {tr  Ml  <_  6} 

Where(H)is  given  by  (1). 

Thus  we  get 

Pr  {tr  S^1  11  <  el  =  6(e)  -  1  [h^A"1)]  c'(e) 

+  0(n‘2) 

on  further  assumption  of  (2)  and  (3)  we  get 

Pr  {tr  S  1  ZZ  <  e}  =  G(e)  +  i  \  a^m.p) 

j=0  3 

5 

V2j(9-")+i  ^  b3(i".P)  9TO*2j(6-“2) 

+  0(n"2)  (6) 

Where  a.j(m,p),  bk(m,p);  j-0,  ....  4,  k*0,  ....  5,  G(e)  and  Gmp(e,w2) 
are  defined  earlier. 
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5.  Sunmary 

Summarizing  we  state  the  following 

Theorem.  Let  Z  *  {Zp  ....  zm>  be  a  p  x  m  random  matrix  where  z^'s 
are  as  defined  In  the  formulation  and  let  nSp  be  a  p  x  P  matrix  which 
follows  a  central  Wlshart  distribution  W(n,p,A"^)  independently  of  Z,  then 
let 

(i )  B’1  A  *  I  +  F  and  |ch1  (F)|  <  1 ,  1-1 ,  p 

(11)  terms  involving  f^f^  are  negligible  where  f^  Is  the  (1 ,j)th  element 
of  F. 

Then  an  asymptotic  expansion  for  percentile  and  c.d.f.  of  T  «  Sn_1  ZZ 
are  given  by  (5)  and  (6)  respectively. 
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6.  Remarks 

(a)  Putting  F  *  0  In  (6)  we  get  the  expression  due  to  Slotanl  [6]  up 
to  the  Indicated  order.  We  note  that  bj(m,p)  terms  vanishes  in  this  case. 

(b)  Putting  M  =  0  In  our  model  we  get  the  expression  given  In  [2]. 

This  should  be  Immediate  If  we  put  M  *  0  and  note  that  In  this  case  J 
in  (2)  reduces  to  corresponding  expression  in  [2], 
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